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ABSTRACT 

I propose a method to detect planets around compact binaries that are strong sources of gravitational radiation. 
This approach is to measure gravitational-wave phase modulations induced by the planets, and its prospect is 
studied with a Fisher matrix analysis. I find that, using the Laser Interferometer Space Antenna (LISA), planets 
can be searched for around ~ 3000 Galactic double white dwarfs with detection limit > 4Mj (Mj ~ 2 x lO'^'g: 
the Jupiter mass). With its follow-on missions, planets with mass > IMj might be detected around double 
neutron stars even at cosmological distances z ^ L In this manner, gravitational wave observation has potential 
to make interesting contributions to extra-solar planetary science. 
Subject headings: gravitational waves — binaries: close — planetary systems 



1. INTRODUCTION 

Since gravitational wave (GW) detectors have omni- 
directional sensitivity, many sources can be simultaneously 
observed without adjusting detectors for individual ones. 
While this might look advantageous for astrophysical stud- 
ies, it also has downsides. Depending on the number of GW 
sources, overlaps of signals in data streams of detectors be- 
come a significant problem, especially in the low-frequency 
regime probed by space GW detectors. For example, LISA 
(Bender et al. 1998) will detect 3000 double white dwarf 
binaries above ~ 3mHz (see e.g. Nelemans 2006, Ruiter et al. 
2007). Without removing the foreground GWs made by these 
numerous binaries, it might be difficult to observe weak in- 
teresting signals, such as extreme-mass-ratio-inspiral (EMRI) 
events. In this respect, extensive efforts are being paid to nu- 
merically demonstrate how well both strong and weak signals 
can be analyzed, using mock LISA data (Arnaud et al. 2007). 

In this paper, I study a method to search for plan- 
ets (more generally sub-stellar companions) orbiting around 
ultra-compact binaries. The proposed approach is to observe 
binaries' wobble motions caused by the planets and imprinted 
as phase modulations of GW from the binaries. This ap- 
proach is close to the eclipse timing method (see e.g. Deeg 
et al. 2008) to detect planets around binaries, and the under- 
lying technique is similar to the planet search around pulsars 
with radio telescopes (Wolszczan & Frail 1992, see also Dhu- 
randhar & Vecchio 2001). As the expected modulations due 
to planets are small, the ongoing numerical efforts for LISA 
have direct relevance to the prospects of the detection method 
proposed in this paper. 

Here, I briefly discuss the significance of this method on 
extra-solar planetary science. In the last 15 years, its rapid 
progress has largely been led by theoretically unanticipated 
discoveries, such as those of the hot Jupiters (Mayor & 
Queloz 1995) or the pulsar planets (Wolszczan & Frail 1992). 
However, at present, observational studies for circum-binary 
planets are in a very preliminary stage (Udry et al. 2002, 
Muterspaugh et al. 2007, Deeg et al. 2008). In addition, im- 
pacts of stellar evolution processes including giant star phases 
or supernova explosions are still highly uncertain (see e.g. 
Villaver & Livio 2007, Silvotti et al. 2007 for recent studies). 
Since ultra-compact binaries such as double white dwarfs are 
end products of stellar evolution, the proposed method to 
search for planets around them would provide us with im- 



portant clues to these unclear issues. While the probability 
of finding a planet around a compact binary is uncertain, the 
large numbers of available binaries (e.g. with LISA) are ad- 
vantageous for various statistical analyses, such as estimation 
of mass distribution of planets by separating information of 
orbital inclination sin/. 

2. PHASE MODULATION BY A PLANET 

To begin with, I discuss GWs from a detached double white 
dwarf binary on a circular orbit without a planet (Takahashi & 
Seto 2002). I write its almost monochromatic waves around 
frequency fg^. as 

ho(t)=Acos[2Trfg,,t+TTfg„f+ipo+DEit)]=Acos[ip(t)], (1) 

where the term (cx fga.) represents the intrinsic frequency evo- 
lution with fgwTobs ^ fgw {Tabs'- observational time < lOyr). 
The term Z)£(f) = 27r/g„,i?£C~' sin0iCos[(/)(f)- (/>j] represents 
the Doppler phase modulations due to revolution of a detec- 
tor around the Sun {Re =1AU) with its orbital phase 0(f) = 
2TT{t + const. The angular parameters iOs,<j>s) are the di- 
rection of the binary on the sky in the ecliptic coordinate. In 
eq.(l) I have neglected amplitude modulation by rotation of 
the detector To determine the direction of the binary, this ef- 
fect is less important than the Doppler modulation DE(t) at 
fgK ^ c/Re ~ ImHz (Takahashi & Seto 2002). In relation 
to this, I do not explicitly deal with the orientation parame- 
ters of binaries. This is just for simplicity. These parameters 
determine the polarization states of the waves. 
The orientation-averaged amplitude of the waves is given 

as 
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with the chirp mass M = m\' M^l {Mi +M2)"'/^ {Mi and M2'. 
two masses of the binary). In this equation, I put the chirp 
mass at = 0.45Mq (Farmer & Phinney 2003) and used the 
distance to the Galactic center r = 8.5kpc as the typical dis- 
tance to Galactic binaries. The matched filtering technique is 
an advantageous method for GW observation and the signal- 
to-noise ratio of the binary is evaluated in the standard manner 
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with LISA detector noise level hj that is within 15% around 
1 .2 X 10-2" Hz-'/^ in the frequency regime 3mHz-10mHz rel- 
evant for the present analysis (Bender et al. 1998). Here I 
assumed that LISA has two independent data streams with 
identical noise spectra. 

When the binary has a circum-binary planet with mass Mp 
and orbital frequency fp, the observed waveform hnjit) has 
an additional phase shift Dp{t) due to the binary's wobble 
induced by the planet, and 1 put the waveform by /im(0 = 
Acos[ip(t) + Dp{t)]. For a planet on a circular orbit, the phase 
shift is given by Dp{t) = pCOSLpp(t) with the orbital phase 
(Pp(t) = lirfpt + ifco phase constant) and the amplitude 
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tal mass of the binary, inclination of the planet's orbit) or 
explicitly 



(6) 



(7) 



^3mHz J V0.33yr-i / ■ 
For a system at a cosmological distance with redshift z, the 
amplitude ^'p is given by multiplying a factor (1 +z) to eq.(|5]l 
with the intrinsic (not redshifted) orbital frequency fp and the 
observed (redshifted) GW frequency fgu,. Note that the red- 
shift z can be estimated from the observed luminosity distance 
(Schutz 1986). As 1 want to know the smallest mass MpSini 
detectable with GW observation and it is easier to find a planet 
with a larger amplitude ^fp, 1 hereafter assume 'i'p-^i 1. Then 
the modulated signal h^it) is expressed as 

hM(t) = ho(t) + hp+(t) + hp. 
with two new components 

hp±(t)=A^p(sm[ip(t)±ippm/2 

= -A^'p(sin[27r(/^„ ± fp)t + 7r/j„f ^ + ip^± (^po])/2. 

A simple interpretation can be made for eq.®. In addition to 
the original signal ho{t) given in eq.(l), motion of the planet 
produces two replicas hp± (smaller by a factor of 5',,/2 than 
hoit)) 3.1 nearby frequencies fg^. ±fp^ fp. Because of the 
coupling with the binary's rotation, the orbital frequency fp 
of the planet is now up-converted into a band that might be 
observed with GW detectors. Here, it is important to note 
that the gravitational wave signal of each replica hp+ or hp- 
itself is described with a nearly monochromatic waveform for 
a standard Galactic binary (including dependencies on angu- 
lar parameters). This fact is important for data analysis, as 
seen later In this paper 1 only study a planet on a circular 
orbit, but this analysis can be straightforwardly extended for 
multiple planets or eccentric orbits that produce other small 
replicas at frequencies /j,„. ± nfp (n = 2, 3, • • •) not only with 
n = 1 (Dhurandhar & Vecchio 2001). 

Based on the simple interpretation of the modulated signal 
hmit), 1 can naively define the signal-to-noise ratio for the two 
small replicas hc± by 
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as for the original one ho given in eq.(3). If the parame- 
ters ao = iA,fgH.,fg^^-,Lpo,0s,4>s) are well determined with the 
strong original one ho, they can be used to estimate the three 
additional parameters a.^ = (/p, ^'p, V'eo) for the small repli- 
cas hp±. The expected observational errors for the three new 
parameters cx^ are evaluated by a 3 x 3 Fisher matrix (see e.g. 
Takahashi & Seto 2002), and I obtain the asymptotic results 
at Tohjp > 1 as 

fA^p 



(9) 



where the suffix "3" represents fitting only three new param- 
eters. 

The actual observational situation is more complicated. For 
example, the frequency resolution is given by ^ T^"^' , and with 
a short observational period Tabs, there must be a significant 
interference between the weak signals hp±(t) and the strong 
one hoit). Furthermore, the triplet {ho,hp±) is need to be 
identified in the presence of thousands of other binaries. To 
study the interference within the triplet, 1 firstly discuss signal 
analysis only with a single binary-planet system and detec- 
tor noises. I numerically evaluated the observational errors 
expected for simultaneous fitting of all the nine parameters 
ao and listed above, and obtained the magnitudes of the 
errors (A^p/^p)^ and (Afp)s for various sets of input pa- 
rameters (suffix "S": simultaneous fitting). 1 found that, for 
a given observational time T^bs ^ lyr, these results depend 
strongly on the orbital frequency f,, weakly on the phase (pco, 
and negligibly on other parameters. As shown in figure 1, 
the errors (A^'p/^'p)^ and (A/p)s become much larger than 

the previous simple estimations (A^p/^p)^ and (Af,)^ for 
frequencies 

fpTobs<2 or |/p-l|7;,ft,<l. (10) 

In the latter regime, two phase modulations Z)£(f ) and Dp{t) 
become highly degenerated. Outside these two bands, the 
replicas hp± are well separated from the original one ho in 
the frequency space, and the simple estimation in eq.(|9l) be- 
comes reliable. In these preferable frequency regimes, the 
mass Mp sin/ can be estimated within 10% error (at the same 
time, the naive SNR Xp > \Q) for planet with Mp sin/ > 
5.4My(//3mHz)-5/3(/p/0.33yr-^)2/3. Here 1 used eqs.©® 
and (|9]), and the following typical parameters: r = 8.5kpc, 
Ml = M2 = 0.52Mq, Tob, = lOyr and hf = 1 .2 x 10-2"Hz-'/2. 

Now 1 study circum-binary planet searches among grav- 
itational waves from other binaries. For simplicity, 1 pick 
up a binary-planet system at / > 3mHz and outside the in- 
terfering frequency regimes (fTOl l. 1 consider the follow- 
ing two-steps data analysis; (i) detecting the individual sig- 
nals ho, hp+ and hp-, and (ii) identifying a triplet com- 
bination caused by a planet. The frequency distribution 
of Galactic white dwarf binaries is modeled as dN /df = 
0.08(A?B/3000)(/^v,/3mHz)-ii/3 [yr-'] with the total number 
Nb ^ 3000 at fg^, > 3mHz (see e.g. Bender et al. 1998). A 
similar density is expected for Galactic AM CVn stars (Nele- 
mans 2006). For observational period Tabs ~ lOyr, the occu- 
pation number 7^^"^^ dN /df of binaries per frequency bin will 
be much smaller than 1 . For a planet search, it is crucial to 
detect replicas hp± whose signals are weak but individually 
fitted with standard Galactic binary waveforms. Identification 
of weak binary signals is currently one of the most important 
topics on LISA data analysis. While the situation is some- 
what different, Crowder & Cornish (2007) demonstrated that 
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Fig. 1 . — Planet search sensitivity around wliite dwaif binaries witli LISA. 
Estimated observational eiTors are presented for the orbital frequency of the 
planet ( (A fp)s'- long-dashed curves) and for the amplitude of the GW phase 
modulation ((A^E'p/'I'p)^: solid curves) induced by the planet. These errors 
are normalized by their asymptotic values ((A^p/>I',,)j ,(A/p)3) derived 
with a simple interpretation for the signal modulation (see eq.(9)). Thick 
curves ai'e for integration period Tofc, = 3yr and thin ones for Tabs = lOyr. It 
is difficult to find a planet with a low orbital frequency fpTots ^ 2 due to the 
poor frequency resolution. Two phase modulations DE(t) and Dpit) induced 
by motions of LISA and the planet degenerate at \(fp - \ )\Tgi,^ <^ 1. Outside 
these two bands the new signal hp±(t) by the planet can be well separated 
from the strong original one /jqCO, and the planet search works efficiently. 
These results depend very weakly or negligibly on source parameters other 
than . 



many (but not all) binaries can be detected down to SNR ^ 7 
(corresponding to Xp/V2 ^ 7 for each replica) even under 
a more crowded condition i.e. a larger occupation number 
Tohs dN/df (see their §4.2 and §4.3). They also showed that 
the Fisher matrix analysis provides a reasonable prediction for 
parameter estimation errors. These results are very encourag- 
ing for a planet search that might reversely provide another 
motivation for ongoing activities for LISA data analysis. 

Next I discuss an outline for identifying a triplet signal by a 
binary-planet system. The first task is to search for a potential 
pair h()-hp+ from a list of resolved binaries, using the fact that 
the pair should have same direction (and orientation) param- 
eters with similar frequencies. The second task is to confirm 
the existence of another replica hp- whose parameters can be 
estimated only with the hQ-hp+ pair. Considering the expected 
binary density dN/df, this discrimination method will work 
well. In this manner the triplet can be identified among other 
binaries with a small extension of the standard Galactic bi- 
nary search. Then coherent analysis can be performed for the 
modulated signal h„, to improve the quality of parameter esti- 
mation. 

For unambiguous detections of planets, other effects that 
produce similar waveforms should be closely examined. 
From the arguments about the triplet structure, it is expected 
that the phase modulation Dp{t) can be easily separated from 
other small modulations at higher frequencies ^ fp that also 



Fig. 2.— The ranges of detectable planets for LISA and DECIGO/BBO 
for typical sets of pai'ameters. In the shaded regions, the mass Mcsini can 
be estimated within 10% error. Around lAU (coiTesponding to fp = lyr"'), 
performance of LISA is degraded due to the degeneracy of two phase shifts 
induced by orbital motions of planet and LISA. 



generate small replicas but with larger frequency differences. 
Meanwhile, because of geometrical nature of gravitational 
wave generation, an observed waveform depends on angular 
parameters describing configuration of a binary. For example, 
it is shown that, for an eccentric binary in the LISA band, an 
triplet waveform can be produced by the periastron advance 
with a frequency difference (9(lyr"') (Seto 2001, Willems et 
al. 2007). But the triplet structure is different from the planet 
case. Precession of orbital plane of a binary (by the spin-orbit 
coupling) can also generate a triplet waveform, and might be 
important for double neutron stars with BBO/DECIGO. But it 
has different amplitude patterns (or equivalently polarization 
states), and has a larger frequency difference. 

In figure 2, I plot the detectable planet on the semi- 
major axis-mass plane. Here I used the relation a = 
\{fp/\yr-^)-'^l^{MT /\Mq)^I^K\] for the orbital frequency fp 
and the semimajor axis a. The planets around /p = lyr"' (cor- 
responding to I.OIAU for Mj = I.O4M0) are excluded due to 
the degeneracy discussed before. 

3. DISCUSSIONS 

The follow-on missions to LISA, such as the Big Bang Ob- 
server (BBO) (Phinney 2003) or the Decihertz Interferometer 
Gravitational Wave Observatory (DECIGO) (Seto et al. 2001, 
Kawamura et al. 2006) were proposed primarily to detect 
stochastic GW background from inflation in the band /„„„ < 
/ ^ fmax with /,„,„ ^ 0.2Hz and ^ IHz. At the lower 
frequency regime / < /,„,„, the foreground GWs by extra- 
Galactic white dwarf binaries would fundamentally limit sen- 
sitivity for GW observation (Farmer & Phinney 2003). In 
contrast, at / > a deep window of GW is expected 
to be opened. To this end, it is crucial to resolve and re- 
move foreground GWs generated by cosmological double 
neutron star binaries (NS-i-NSs) whose estimated merger rate 
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is ^ 3 X lO^yr"'. In addition to NS+NSs, there might be 
double black hole binaries or black hole-neutron star bina- 
ries, while their merger rates are highly uncertain. Here I 
provide a brief sketch for planet search around cosmologi- 
cal NSh-NSs with the follow-on missions. I fix masses of 
NSh-NSs at Ml = M2 = 1.4Mq. 

In the observational band [/„!,■„, /max], a NSh-NS is on its 
final stage before merger. The time left before the merger 
is l(/g„./0.2Hz)"^/^(l yr that severely limits the ob- 
servable orbital frequency fp of a planet. Using the restricted 
1.5-order post-Newtonian waveform (Cutler & Harms 2006), 
I evaluated the expected observational errors in the scenario 
that all the parameters are simultaneously fitted, including two 
phase shifts DE{t) and Dp(t). For various sets of input pa- 
rameters, I examined the observational error for the amplitude 
^tp i = ^p|/j„=iHz, and found that by observing at least three 
orbital cycles (namely f,, > f„, = 3{f,„i,j0.2f/\l+zf''yv-') 
the relative error is given as 
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with the signal-to-noise ratio SNRo for the observed NSh-NS 



Here I assumed a nearly flat noise spectrum (in units of 
Hz"'/^) in the band [/„„■„, /mat] (Phinney 2003, Kawamura et 
al. 2006). For a given orbital frequency /p and signal-to-noise 
ratio SNRo, the mass resolution is better than the previous 
results for LISA. This is because of the higher frequencies 
fgw used in the present case. For fp < /,/, (less than three 
orbital cycles in the observational band), the observational 
error A{M p sin i) /{Mp sin i) becomes significantly larger than 

Due to a limitation of estimated computational power avail- 
able at the time of the follow-on missions ^^2025, the min- 
imal noise level of detectors required to remove NSh-NSs 
corresponds to SNRo ^ 100 for NS+NSs at z = 1 (Cutler & 
Harms 2006). For z = 1 the critical orbital frequency becomes 
/,;, = 19yr"' (semimajor axis ~ 0.23AU for Mr = 2.8Mq), 
and the mass Mc sin i can be measured within 10% error for 
a planet with MpSini > \.2(fp/\9yx~^fl^{SNRQ/l00y^Mj. 
The range of detectable planets is shown in figure 2. If de- 
tected at z ^ 1, the planet is ^ lO*" times as distant as those 
currently found in our galaxy. Note that the estimated merger 
rate of NS+NSs around z ^ 1 is ~ lO^yr"'. The bottom edge 
of the shaded region moves to (0.39AU, Q.52Mj) for NS+NSs 
atz = 0.5. 

I would like to thank an anonymous referee for helpful com- 
ments to improve the draft. 
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